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Abstract 

The dynamics of a "peeling front" or an elastic line is studied under creep (constant load) condi- 
tions. Our experiments show an exponential dependence of the creep velocity on the inverse force 
(mass) applied. In particular, the dynamical correlations of the avalanche activity are discussed 
here. We compare various avalanche statistics to those of a line depinning model with non-local 
elasticity, and study various measures of the experimental avalanche-avalanche and temporal cor- 
relations such as the autocorrelation function of the released energy and aftershock activity. From 
all these we conclude, that internal avalanche dynamics seems to follow "line depinning" -like be- 
havior, in rough agreement with the depinning model. Meanwhile, the correlations reveal subtle 
complications not implied by depinning theory. Moreover, we also show how these results can be 
understood from a geophysical point of view. 

PACS numbers: 
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I. INTRODUCTION 



Creep is one of the fascinating topics in fracture for a physicist: the deformation and final 
fracture of a sample follow empirical laws with a rich phenomenology. It is expected that 
there are similarities and differences with "static" fracture encountered in brittle materials 
such that so-called "time-dependent rheology" is not relevant [1]. However, the phenomenon 
of creep is visible in most any setting regardless of whatever a tensile test might indicate 
about the typical material response. A particular scenario where one can study creep is 
the advancement of a single crack under a constant driving force. One can study this in 
simple paper sheets, and for quite some time it has been noticed that this involves statisti- 
cal phenomena, an intermittent response which could be characterized by " avalanches" , in 
particular of acoustic emission (AE) events 0, ^l, Q, 5]. 

A particular experiment we analyze in this work is related to the dynamics of a crack 
line as it moves through a sample, largely constrained on a plane. This can be achieved in 
the case of paper in the so-called Peel-In-Nip (PIN) geometry (see below for a description). 
The tensile case has been already reported in Ref. p and an early account of the creep 
results published as Ref. The mathematical description of the line is a crack position 
h(x,t), where h is the position coordinate along the direction of line propagation and x is 
the coordinate perpendicular to h. On the average, the crack moves with the creep velocity 
v (h = vi). 



The problem has here, as in other such examples (the Oslo plexiglass experiment [9l. llOj). 
three important ingredients: randomness in that the peeling line experiences a disordered 
environment coming from the fiber network structure, a driving force K e ff or a stress in- 
tensity factor, and the self-coupling of the interfacial profile h. In this particular problem, 
it takes place via a long-range elastic kernel 11], expected to scale as 1/x or as k in Fourier 
space. 

For a constant force K e ff the dynamics exhibits a depinning transition, of non-equilibrium 
statistical mechanics. This implies a phase diagram for v(K e ff). The crack begins to move 
(v > 0) at a critical value K c of K e ff such that for K e ff > K c . In the proximity of K c 
the line geo metry is a self-affine fractal with a roughness exponent (. The planar crack 



problem 
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13| has been studied theoretically via renormalization group calculations and 



numerical simulations, and via other experiments as noted above. The roughness exponent 
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of theory Ctheory ~ 0.39 has traditionally been considered to be absent from experiments 



id, 



141 ] , but recent results of Santucci et al. imply that the regime might be visible upon 
coarse-graining. Imaging experiments prove in that case that as expected the line moves in 
avalanches, and the avalanche size s distribution seems to have the form P(s) ~ s -i-6—-i.7 

Here we look at the scenario of creep for the PIN geometry. This subject is such that 
ordinary "fracture creep" and the particular scenario related to depinning transitions coin- 
cide. The creep of elastic lines becomes important for K e ff < K c since thermally assisted 



movement due to fluctuations takes place with a non-zero temperature 15|, [l6|, In usual 



depinning, it is assumed that thermal fluctuations nucleate "avalanches' which derive their 
properties from zero-temperature depinning, and the avalanches then translate into a finite 
velocity v creep > 0. There are two interesting differences in the fracture line creep to other 
such in depinning. First, the line elasticity is non-local, and second, in materials (such as 
paper here) where there is no healing, the line motion is irreversible, there are no fluctuations 
in metastable states as in the case of magnetic domain walls, for instance. 

In this scenario, the creep velocity becomes a function of the applied stress intensity factor 
and the temperature, y rrPPp = v creep (K e ff,T). As creep takes place via nucleation events 
over energy barriers [la ], the description of those barriers is of fundamental importance. One 
can show by scaling arguments and more refined renormalization group treatments that the 
outcome has the form of the following creep formula 

Vcreep ~ exp ( -C/K£ ff ) . ( 1 ) 

This gives the relation to the driving force K e ff using the creep exponent, \l. The value of 
the exponent depends on the elastic interactions and the dimension of the moving object (a 
line), and we expect 

fi = B/v= — . (2) 

a — Q 

The exponents 6, v, and ( denote the energy fluctuation, correlation length, and equilibrium 
roughness exponents. All these exponents are functions of a, the fc-space decay exponent of 
the elastic kernel. For long range elasticity, one would assume a — 1. 

The fundamental formula of Eq. ([2]) has been confirmed in the particular case of 1+1- 

n □ 

dimensional domain walls and other experiments [18|, [19J. We have ourselves reported on 
results, which show an inverse exponential dependence of v CTeep {m) ~ exp (— 1/m), where 



m is the applied mass in the experiment (see below), as is appropriate for non-local line 
elasticity with an equilibrium roughness exponent of ( = 1/3. In the current work we go 
further by two important steps. First, we consider creep simulations of an appropriate non- 
local line model and compare the avalanche statistics and v(m) to those from the experiments 
(see Fig. [T]for an example of the activity timeseries from an experiment and a simulation). 
Then, we ask the fundamental question: what can be stated of the correlations? This relates 
to the timeseries of released energy, to aftershock rates and we present extensive evidence. 
The experimental signatures show subtle correlations that are rather different from what one 
would expect from the (depinning) creep problem with non-existing avalanche to avalanche 
correlations. 




80000 85000 90000 95000 

t 



FIG. 1: Activity as a function of time inside a given time window (a) for the creep experiment 
with 410g load, and (b) for simulations with / = 1.87 and T p = 0.002. In both cases we neglect 
the duration of the avalanche and we only take into acount the starting time and the size of each 
avalanche, obtaining a data series {ti, E{\ for the experiments and {ti, Si} for simulations (definition 
of Si is given in Section III Bp 

The structure of the rest of the paper is as follows. In the next Section, we discuss 
the experimental setup and the simulation model. Section III shows results on v (m) both 
from experiment and simulation. In Section IV we present data on avalanche statistics 
again comparing the two cases. Section V offers an extensive analysis of correlations by 
using a number of techniques to look at the experiment. Finally, Section VI finishes with 
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Conclusions and a Discussion. 



II. METHODS 
A. Experiment 

In Figure [2] we show the apparatus [6]. The failure line can be located along the ridge, in 
the center of the the Y-shaped construction formed by the unpeeled part of the sheet (below) 
and the two parts separated by the advancing line. Diagnostics consist of an Omron Z4D-F04 
laser distance sensor for the displacement, and a standard plate-like piezoelectric sensor Q]. 
It is attached to the setup inside one of the rolls visible in Fig. [21 and the signal is filtered and 



FIG. 2: Experimental setup for peeling experiment. The paper (white) is peeled between two 
cylinders (copper) separated by a few millimeters. The driving force is generated by a larger 
hanging weight (black). A smaller weight adjusts the peeling angle. The AE and distance data are 
collected by piezo transducer (red) and a laser sensor (gray). 

amplified using standard techniques. The data aquisition card gives us four channels at 312.5 
kHz per channel. We finally threshold the AE data. The displacement data is as expected 
highly correlated with the corresponding AE, but the latter turns out to include much less 
noise and thus convenient to study. For paper, we use perfectly standard copy paper, with 
an areal mass or basis weight of 80 g/m 2 . Industrial paper has two principal directions, 
called the "Cross" and "Machine" Directions (CD/MD). The deformation characteristics 

n 

are much more ductile in CD than in MD, but the fracture stress is higher in MD [20]. We 
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tested a number of samples for both directions, with strips of width 30 mm. The weight used 
for the creep ranges from 380 g to 450 g for CD case and from 450 g to 533 g for MD case. 
The mechanical (and creep) properties of paper depend on the temperature and humidity. 
In our setup both remain at constant levels during experiments, and the typical pair values 
for environment is 40 rH and 26 °C. 



B. Simulations 



We want to simulate the evolution of a discrete long-range elastic line of size L in a 
disordered media. The line is characterized by a vector of integer heights {hi, . . . hi) with 
periodic boundary conditions, although the experiment does not present periodic bound 
conditions. 

The long-range elastic force [2JJ acting on a string element is given by 

7T\ 2 \~\ hi — hi 



f elastic = k M V "J " (3) 

\LJ ^ S in(^^7r) 

where all forces on all sites can be computed in a L log L operations using a fast-Fourier- 



transform (FFT) algorithm [22|. Simulations are done using k = 0.01 and L = 1024. 
The random force due to the quenched disorder may be obtained from a standard normal 
distribution, i.e a Gaussian distribution with zero mean and a variance of one f^ andom = 
N(0, 1). Then, the total force acting in a given element of the string is fi = ff astlc + 
^random _|_ j ^ w h. e re / is the external applied force. 

At this point, we need to introduce a dynamics which mimic the experiment evolution. 
A basic characteristic of the experiment is that it is completely irreversible, so the dynamics 
has to include this important feature. We consider a discrete time evolution and the discrete 
dynamical rule [3] is given by: 

Kit + 1) - Kit) = Vi (t) = 9\fi] t = 0, 1, 2 . . . (4) 

where 9 is the Heaviside step function. Then we apply the following procedure: 

1. Start at t = with a flat line located at h = setting K = Vz. 

2. Compute the local force (fi) at each site and using the dynamical rule (Eq. @| compute 
the local velocity of each site. We can define the velocity of the string for this time, 
as v(t) = xEf=i 



3. Advance the sites according their local velocities 

4. Generate new random forces for those sites that have been advanced. 

5. Go to step (2) and advance the simulation time by one unit. 

This evolution shows a depinning transition at f c ~ 1.88 in which the velocity of the line 
v(t — > oo) > when f > f c and v(t — > oo) = when / < f c . 

In order to simulate the creep evolution of the string we use an external force below the 
depinning threshold, and when the line gets stuck we let thermal fluctuations play a role. 
We scan all the sites and set Vi{t) = 1 with a probability p = exp an d Vi(t) = with 
a probability 1 — p, where T p is proportional to temperature. This can trigger an avalanche 
which will have a finite duration T since the system is below the depinning treshold. We 
define the avalanche size as S = Y1t v (^)- ^ we consider small enough temperatures com- 
pared to the typical internal forces, the avalanche needs some time to be triggered, which is 
defined as the waiting time r. We define this waiting time as the time between the end of 
an avalanche and the starting time of the next one. 




FIG. 3: Velocity of the long-range elastic string as a function of simulation time (dotted line). The 
vertical and solid lines represents the signal S(t) ploted in figure [TJ(b) . Avalanche properties are 
also shown: is the waiting time, T{ is the avalanche duration duration, and Si is the avalanche 
size. 

In summary, this long-range elastic line model in the creep regime has an avalanche- like 
behaviour. Each avalanche is characterized by three quantities: Waiting time r, duration 

7 



T, and size S (see Fig. [3]). Moreover, we observe that for long times, when the steady state 
is reached, durations are small compared to waiting times, for that reason we can simplify 
the signal just taking into account the starting time of the avalanche and its size. 



III. CREEP VELOCITY 



The main data about both simulations and experiment on the creep velocity are shown 
in Figure HI The prediction of Eq. (TjQ) is that the velocity is exponential in the effective 
driving force. In the case of the experiments at hand, we face the problem that we do not 
know the average fracture toughness (K) empirically. It depends on the loading geometry, 
and on the material at hand. There are estimates for similar papers in the literature in the 
mode I case (see e.g. Ref. 



231 ]) which indicate that the value of (K) is lower by at least 
a factor of two compared with the actually used loads. One can try to work around the 
problem by guessing (K) and checking how that affects the apparent functional relationship 
of v vs. the reduced mass m e ff = m — (K). In the range of physically sensible values of (K) 
the velocity exponential behaviour does not change and thus we can take m e ff = m . 
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FIG. 4: The creep velocity vs. the inverse of the applied force or mass, m e ff = m. Inset: creep 
velocity vs. / for the simulation model for two different temperatures. 

From the figure we may conclude that the effective creep exponent \x ~ 1, though there 
is variability among the data sets. One of the data sets (black circles) shows some slight 
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curvature. The main finding, interpreted via Eq. (|2j) then indicates that the effective 
roughness exponent C ~ 1/3, which is the expected equilibrium value for a long-range elastic 
problem with a = 1 7]. 

The numerical simulation data agree qualitatively with the exponential d ecay except very 



close to the depinning transition. According to the creep formula (see e.g. 17]). we should 
expect that the velocity of the long-range string was 



v(f,T p 



exp 



(5) 



However, it appears that slope as a function of the temperature is not exactly the expected 
one. One reason is that the model is simplified: we only let thermal fluctuations act when 
the string gets stuck so avalanche nucleation during an avalanche is neglected. This may be 
of importance very close to f c and for long avalanches. 




FIG. 5: Histogram of a normalized velocity obtained from discretized distance data. Velocity, v, 
is an average in a 0.5 s time window. < v > is an average over experiments with same weights. 

The exponential average creep velocity can most directly be compared with the measured 
velocities from the distance sensor over short time-spans. Figure [5] shows the probability 
distributions P{v) for a very large number of different experiments, for the v = Ah/ At with 
At = 0.5 s. The general trend shows clear stick-slip characteristics in the sense that the local 
velocities vary with a power-law -like fashion. The typical slope of the data is about -2.3 
though a more detailed look indicates that there is a tendency for the exponent to change 
with m and with At (increasing both decreases the slope). It is an interesting question 
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of how this locally time-averaged velocity is related to the average creep velocity, and the 
avalanches that contribute to it, somewhat hindered by the relative large fluctuations in the 
distance sensor - for which reason we resort in the detailed avalanche dynamics studies to 
the AE. 

IV. STATISTICAL DISTRIBUTIONS 

Next we consider the statistics of the AE timeseries from the experiments as signatures 
of the intermittent avalanche activity in the system during creep. In our setup, we face 
the problem that direct imaging of the front dynamics is if not impossible then difficult to 
realize. Thus we take the AE data up to be scrutinized as detailed information. It can be 
studied from the viewpoint of the correlations of the creep avalanche activity but the finer 
details thereof are left to the next Section. Here, we consider the averaged distributions of 
three quantities: i) AE energy E in experiments and avalanche size in simulations assuming 
that the avalanche energy is proportional to its size E ~ S ii) waiting times r both for 
experiments and simulations and iii) avalanche durations T only for simulations, because 
experimental avalanches have very short durations and can be neglected. 

, „0 
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FIG. 6: Energy distributions for the tensile experiment (circle), for the creep experiment (square), 
and for the simulations (triangle up). For the simulations we are ploting the histogram of the 
avalanches sizes {Si}. We can consider that the energy of an avalanche is proportial to its size, so 
S t ~ Ei}. 
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FIG. 7: (a) Waiting time distributions for the tensile experiment (circle), for creep experiment 
with 410g (squares), and for the simulations with / = 1.87 and T = 0.0020 (triangles up), (b) 
Normalized waiting times for different creep experiments. 

Figure [6] shows three cases of the avalanche size distributions. We compare the creep data 
for one mass m to a similar dataset for a tensile experiment done at a constant average front 
velocity 6j. Moreover data is included from the creep model for the parameters shown in the 
caption. The normalization of the data for the experiments is such that the E min has been 
scaled to unity. Recall that the events are restricted in size from below by a thresholding 
applied to the original AE amplitude signal A(t), from which the events are reconstructed. 
We can observe that the effective power-law exponents of the experimental data are ~ 1.6 
for the creep and ~ 1.8 for the tensile cases, respectively. These are very close to each other, 
while the simulation data results imply ~ 1.4 not very far from the experimental values. 
We also can observe that there is no evident cutoff in any of them (the bending in the 
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simulations case is a finite size effect). These data can be compared with the Oslo plexiglass 
experiment where for the avalanche size distribution the value of (3 = 1.6 ± 0.1 has been 
found QQ. 




T 



FIG. 8: Avalanche duration distributions for Tp = 0.002 and three different forces. For the case 
with / = 1.89 we are above the depinning treshold. 

The waiting times are reported in Fig. [7J For all the three cases -P(t) is broad. In the 
tensile case, it is known that there appears to be a "bump" in the distribution, or a typical 
timescale. This is absent from the creep one. It is interesting to note that here the simulation 
model agrees rather well with the creep case. For larger m it is possible that the waiting 
times start to look more like the tensile case. We also present the scaled distributions for 
all the experiments. Later, in the next section, we discuss the attempt to link this to a 
background plus correlated, triggered activity. 

Finally, in Fig. [8] we show the avalanche durations from the simulations. In the case of 
the experiment this is more complicated due to the fact that the actual amplitude signal 
is convoluted via the preprocessing electronics and the response function of the piezos with 
which the AE is measured. Later we present some examples of the outcome, but here we 
just discuss the clear-cut case of the simulations also since they give an idea about what one 
might see in the experiment, ideally. The main points that one learns from the Figure are 
that a true power-law-like P{T) ensues only at the proximity of the f c . For values higher 
or lower than that the shape of the distribution changes, in particular such that not only a 
cut-off appears but also the clear power-law character starts to disappear. 
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V. MEASURES OF CORRELATED DYNAMICS 



A. Correlations 



Next we look at the detailed temporal structure of the AE signal. The main question is 
whether the creep activity exhibits interesting features that would in particular differ from 
the theoretical expectations - based on elastic line depinning the inter-avalanche correlations 
should be expected to be negligible. 
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FIG. 9: The autocorrelation function of the averaged event energy in paper peeling under creep 
and tensile loading modes. Comparisons to the randomized data are also included. The numerical 
data and the corresponding randomized data are not distinguishable. 

In Figure [9] we show the autocorrelation function R(t) of the event energy time series. 
The autocorrelation function is defined as: 

R{u) - — wrw) — {) 

where E t is the energy of the AE signal at time t and (E) is the average value of the energy. 
Et is defined as a sum of squared amplitudes of the AE signal in the time interval [t, t + At]. 
The length of the interval At is chosen to be 10~ 3 s in the tensile and 10~ 5 s in the creep 
peeling experiment in order to capture the correlations in both cases. 

When compared to paper peeling experiments under a constant strainrate, the correlation 
decays at a much faster rate than in creep peeling experiments. In that case, the existence 
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of a slow decay might be taken to be connected to the fact that there is a typical scale 
in the waiting time distribution which is not the case for creep, seemingly. The functional 
form of the shown case of a logarithmically decreasing autocorrelation function is R(t) = 
—0.3 — 0.08ln(t). The data are also compared to a randomized timeseries, and one can see 
that the correlations disappear. For the simulated data the autocorrelation function shows 
no difference to a randomized signal. All in all these results imply that there are contrary 
to theoretical models temporal correlations, albeit in creep on a very short timescale. 

In Figure [10] we show an envelope event form for different events with different event 
energies. We see an exponential decay for the tail of the event, but the event envelope 
becomes more extended in time when the energy of the event is larger. Typical events extend 
up to 0.5ms, which corresponds to decay of the correlation in the Fig. A correlation up to 
time-scales larger than the typical event length is only seen in the strain-controlled peeling. 

I i i i i i i 
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FIG. 10: Squared amplitude of an event averaged over all events in the creep peeling experiment. 
The average is taken over events with different sizes and the size is indicated as different colors in 
the figure. 

Since the timeseries of AE is so intermittent it is better to concentrate on measures 
that consider directly the avalanches. In Figure [TT] we depict the averaged energy as a 
function of a silent time before the event from paper peeling experiments in creep. The 
event energy is in many datasets slightly correlated to the waiting time before the event. 
This correlation disappears if one considers the opposite case of the waiting time after the 
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event. The suggested interpretation is that the elastic fracture line apparently as a physical 
system ages before a large event, while there is no real dependence of the waiting time on 
the energy dissipated in the previous event. 
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FIG. 11: Averaged energy as a function of silent time before the event with weight 410g. 

The difference in the autocorrelation between the creep and tensile peeling experiments 
might be attributed to the forcing the line to move in the latter, which induces a "fiber- 
scale" to results. This is also supported by observing the waiting time distribution, where 
the pdf deviates from a power-law. 

In paper peeling we study the clustering of events by computing the correlation integral 
C(AT), that is the probability that two events are separated smaller time than AT. The 
correlation integral is given by: 

C(AT) = 2 _ 6 ( AT " (*i " **)] > ( 7 ) 
where N is number of events in the experiment and £j is the event occurrence time. 



Correlation integrals 24j are shown in the Fig. [12] for the peel creep experiment. If the 
probability of the event occurrence is equal for every time interval, then one can assume that 
correlation integral increases as C(AT) ~ AT. We see a power law AT 0,9 in sufficiently 
large times, but when the distance of events approaches the experiment length we see small 
deflection in the curve. At temporal scales of the order of 10 _2 s we see deviation from the 
power law behaviour, which indicates event clustering. 
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FIG. 12: Correlation integrals for the creep in peeling experiment. 
B. Seismicity - cascading occurrences as a model for the experimental data 

In this part we will show how fracture in heterogeneous material, such as line creep 
in paper peeling, behaves, in time, similarly to the rupture at the Earth scale, e.g. the 
earthquakes driven by plate tectonic deformation. 

From seismology it is known that seismicity can be described by two processes: the 
background seismicity and the triggered events. The first one is modelled as a homogeneous 
Poisson process, while the second one as a power law decay of seismic rate following the 



occurrence of any event, e.g. the Omori's law 



2a 



26 



27|: 



R = fi + J2 X i( t ) ( 8 ) 

t<ti 

The first term in the right hand side of Eq. [8] is the background seismicity, while the second 
term is the correlated part of the seismicity, that is, the superposition of time- dependent 
series of triggered seismicity following any event. The triggering process of the latter is 
reproduced by models of cascading effect for earthquake interactions, i.e. ETAS (Epidemic 



Type Aftershock Sequence) model 25|, |26j, |28|]. This stochastic point process is based on 



the Gutenberg-Richter law for energy distribution and Omori's law for time distribution of 
seismicity rate. According to this model, the rate of aftershocks triggered by an earthquake 
occurring at time £j with magnitude Mj is given by: 

= ^0_i OQ(Mi -mc) ( 9 ) 

v ; (c + t-U)p v ; 
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where K , a, c and p are constants and mc is the completeness magnitude of the catalogue. 

The total earthquake rate of Equation (jSJ) is therefore the sum of all preceding earthquakes 
(triggered directly by the background events or indirectly by previous triggered events) and 
the constant background rate p,Q. This model reproduces most of the statistical properties 
of earthquakes, including aftershock and foreshocks distributions in time, space and energy 
27|. 

Figure [T3l illustrates the average acoustic event rate following any event for the peel creep 
experiments (load m = 409 g). It is reminiscent of Omori's law for tectonic seismicity, where 
we can observe the power law decay representing the cascade of aftershocks following an 
event. For times greater than 10~ 2 seconds, the event rate keeps constant, at the background 
rate level. The exponent of the power law decay of event rate is equal to 1.5 ± 0.1. This is 
fairly close to what would fit the experimental P(t) demonstrated in Fig. 
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FIG. 13: Event rate following events in paper peel creep experiments with m = 409g. Time t = 
is the target event occurrence. Aftershock rates are averaged within each magnitude class of target 
event (blue line: 4.2 — 5.72; red line 5.72 — 7.2). We compute the magnitude class M = logio(EM) 
where Em is the energy of the target event. All magnitude classes are averaged together (thick black 
line). Correlation between events is characterized by a power law decay of the activity after the 
target event. The time for which events are correlated is a function of the target event magnitude, 
as well as the number of triggered events (see Equations ([8]) and ©). The observed duration of 
the aftershock sequence is bounded by the level of the background uncorrelated constant rate. 

The AE, triggered by line creep in paper peeling, is characterized by power law distri- 
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bution on energy (Fig. [6]) and power law relaxation of aftershock rate (Fig. [T3|) . ETAS 



style models reproduce 
conditional mainshocks 



:hese macroscopic patterns, including foreshocks as aftershocks of 
. Corral 



29j shows that the inter-event time probability density 
for such kind of ETAS model for event occurrences follows a gamma distribution, according 
to: 

p( T ) = Ct^ 1 exp(-r//3) (10) 

where r is the normalized inter-event time obtained by multiplying the inter-event time St 
with the earthquake rate A, that is r = StX. 

Molchan 30j showed that, in agreement with Equation (flOj) . the distribution decays 
exponentially for large inter-event times and that the value 1/(3 is the fraction of mainshocks 



among all seismic events. According to Hainzl et al. 31[, 1//3 is a regional quantity, allowing 
for non-parametric estimate of the background rate in a specific process. In order to simulate 
the AE properties of the creep fracture experiment (m = 409g), we tuned an ETAS model to 
fit the estimated percentage of background activity of real data. One must notice that robust 
inversion of ETAS model parameters is not yet available. Figure [TH shows the comparison 
between inter-event time distributions of a synthetic catalogue generated by ETAS model. 
Both, simulations and data inter-event time distributions fit a gamma distribution. Other 



possible data fittings are possible 32], but this lies outside our aim of comparison between 
data from paper peeling and ETAS simulations. The fit may underestimate here (see also 
Fig. [7b ) slightly the exponent of the power-law part of the waiting-time distribution. In 
any case, the relevant exponent here is definately smaller than in the case of rock fracture 



m (p = i.4). 

To summarize, line creep in paper peeling at a scale of ~ 10 _1 m and ~ 10 2 s triggers 
brittle creep damage that seems to share the same generic temporal properties than the ones 
observed for tectonic seismicity at scales of ~ 10 6 m, ~ 10 2 years. These properties can be 
reduced to a rough constant seismicity rate with bursts of correlated activity, contemporary 
to power law distribution of event sizes and (short-time) inter-event times. Estimates of 
Omori's law exponent suggest a faster relaxation for the paper peeling case than for Earth 



crust response to tectonic loading, p equal to 1.4 and 1 respectively [26]. The portion of 
uncorrelated events suggests a slightly lower triggered event rate in paper peeling than in 
the Earth crust deformation. Estimations of the background portion of AE did not show 
any sensitive dependence on the applied loading. Whether the difference between paper 
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FIG. 14: Inter-event time probability distribution for experimental dataset (thin red curve) and 
synthetic catalogue generated by ETAS model (thin black curve). Dotted thick curves are gamma 
distribution fits to data and ETAS model (red dotted line for the real data and black dotted 
curve for ETAS). Estimations of background fraction of events according to Hainzl et al. (2006) 
technique are close together (23 — 25%) for both data and simulation. ETAS parameters are: 
p = 1.4, K = 0.09, a = 0.9 (n = 0.9), b = 1, c = 0.001s. 

experiments and earthquakes come from experimental conditions or fracturing mode (i.e. 
tensile, creep or compression) remains an open question. For earthquakes no change in 
relative portions of background and triggered activity is resolved for compression, extensional 
or shear tectonic settings. 

VI. CONCLUSIONS 

We have overviewed a simple creep experiment which uses paper and can be studied to 
investigate planar crack propagation in a disordered medium. The information that one can 
obtain and then compare to relevant theory extends from the average front velocity to details 
of the spatiotemporal dynamics. We have also for a comparison studied a classical non-local 
elastic line model under creep conditions. This shows similar features to the experiment: an 
exponential dependence of the creep velocity on the applied force or mass or stress-intensity 
factor. 

The typical statistical distributions are power-law -like in particular for the event en- 
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ergy/size. It is perhaps useful to recall that the waiting time distribution is quite broad. 
There is currently no understanding as to why, in particular one should note that the current 
experimental setup allows to study this issue in a steady-state unlike in most other fracture- 
related creep tests. In general as such distributions are regarded the line creep model agrees 
at least qualitatively with the experimental data. Our results are also in line with other sim- 
ilar planar crack data (though these are obtained usually in the constant-velocity ensemble, 



not in creep 



AH, 



Looking in more detail at the correlations of the activity, differences transpire however. 
The experimental AE events show subtle correlations via the autocorrelation function, via 
the waiting times before events, and via the Omori's law. All these measure different aspects 
of the avalanche activity, and in all the cases the model differs in its behavior. Here, we lack 
completely theoretical understanding, in particular as regards such a quantitative measure 
as the Omori exponent. It is interesting to note that geophysics -oriented analysis methods 
produce results in agreement with observations from tectonic activity. Here again the steady- 
state character of the experiment at hand is of utility. 

In the future such experiments and such comparisons can be used to study several different 
aspects of avalanching systems, creep fracture, and models for line depinning. A particularly 
pertinent question is for instance whether rate-dependent processes in the material at hand 
modify the kinetics of the creep in some suitable way that still maintains the creep vs. force 
-relation intact. We shall ourselves attempt a more careful study of the creep model, and 
analyze how its correlation patterns could be matched with the experiment. 
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